Constraints on a long-range spin-independent interaction from precision atomic 

physics 
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Constraints on a spin-independent interaction by exchange of a neutral light boson are derived 
from precision data on the electron anomalous magnetic moment and from atomic spectroscopy of 
hydrogen and deuterium atoms. The mass range from 1 eV/c^ to 1 Mev/c^ is studied and the 
effective coupling constant a' is allowed below the level of lO"^'^ — 10~^^ depending on the value of 
the boson mass. The mass range corresponds to the Yukawa radius from 0.0002 nm to 20 nm, which 
covers the distances far above and far below the Bohr radius of the hydrogen atom. 
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I. INTRODUCTION 

Precision tests of quantum electrodynamics for bound 
states and free particles allow us to verify various ad- 
vanced methods of measurements and calculations and 
determine precision values of certain fundamental con- 
stants, such as the Rydberg constant R^o and fine struc- 
ture constant a. The overall consistency of the results 
obtained for those constants [l|, proves, in particular, 
that we can consistently describe various fundamental 
phenomena in a broad range of distances and energies [21 . 

Meanwhile, an introduction of new physics could affect 
various scales with different strength and violate the con- 
sistency mentioned. If such a problem is not observed, 
one sets various constraints on possible new physics. 

In particular, various unification theories suggest ad- 
ditional particles (see, e.g., |3-[a|)- ^ specific kind of 
such particles is a light neutral particle weakly interact- 
ing with conventional matter consisting of electrons and 
nucleons. Indeed, 'lightness' in terms of particle physics 
ranges from below 1 eV/c^ to above 1 GeV/c^. Stable 
neutral particles of this kind can be also considered as a 
candidate for the dark matter Q. 

In a broad range of distances one can consider an addi- 
tional new interaction as a modification of the Coulomb 
interaction such as 



a acff(r) 
r r 



-,-Ar 



(1) 



A similar substitute can be also written in the momentum 
space. 

The particle is not necessarily stable in our considera- 
tion. Generally, one should rather expect that a particle 
coupled to charged particles, decays into a few photons, 
while a particle coupled to leptons, could also decay to 
neutrinos. (Decay into a pair of massive charged leptons 
can be forbidden for light particles because of their light- 
ness.) The substitute ([1]) is valid as long as the decay 
width is much smaller than the mass A. 
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Comparing values of the fine structure constant a 
obtained from experiments with different characteris- 
tic distances and momenta, one can check whether the 
Coulomb coupling constant is really a constant. (We 
have in mind that the vacuum polarization corrections 
responsible for the 'running coupling constant' have been 
already taken into account.) 

The most precise measurements of various physical 
quantities are often aiming at the determination of the 
values of fundamental constants. The related precision 
data open a certain window of opportunities to verify 
the constancy of the Coulomb coupling constant. There 
are basically three important scales, which contribute to 
precision QED-related experiments. 

The scale of distances related to the Bohr radius 
flo = l/(aTOec) is studied within hydrogen spectroscopy 
with transitions involving the ground and the low ex- 
cited states (Is, 2s). The distances of about lO'^ao are 
accessible in experiments on high Rydberg states in the 
hydrogen atom (n ~ 30), while the Compton wave length 
of the electron Ac = h/{mec) is a characteristic distance 
involved in calculations of the electron anomalous mag- 
netic moment. 

Altogether, the related distances vary from 2 x 10^^'^ m 
to 2 X 10^^ m and related masses range from 1 eV/c^ to 
1 MeV/c^. Indeed, in the case of the completely relativis- 
tic calculation of a correction to the electron anomalous 
magnetic moment we cannot rely on a static component 
in ([1]), but should deal with an interaction describing an 
exchange by an intermediate particle. 

Let us now consider precision QED-rclated experi- 
ments for these three areas step by step. 



II. ANCHOR EXPERIMENTS: 

MEASUREMENTS ON LOW STATES IN THE 
HYDROGEN ATOM (DISTANCES ~ (1 - 4)ao) 

A massive set of data on hydrogen and deuterium 
spectra is available for low excited states related to 
physics from one to a few Bohr radii. They are used 
here as an anchor for further comparison with other dis- 



tance/momentum scales. We remind that the Bohr ra- 
dius ag — 0.53 X 10~^° m is a characteristic distance for 
the electron orbit at the ground state. For the 2s state 
the characteristic distance is ^ 4ao. 

Throughout the paper the relativistic units are applied 
in which h = c = 1. In these units ap = l/ame- For the 
Yukawa radius equal to the Bohr radius, the related mass 
of the intermediate particle A is 3.5 keV, while the radius 
of 4ao corresponds to 1 keV. 

To understand the procedure of evaluation, one can 
have in mind a simple picture with the non-relativistic 
Schrodinger equation, because all corrections beyond 
such approximation, which are due to QED and relativis- 
tic effects, are well under control and can be introduced 
when necessary [ll, 0, H ■ 

Spectroscopic data, which are the most important sta- 
tistically, are related to the Is — 2s transition [9|, 
to the 2s — ns/d transitions for n = 8, 10, 12 
remind that theory of hydrogen levels can be expressed 
in terms of the Rydberg constant i?oo and the proton 
charge radius [ll, 0; @1 ■ Those constants are not known 
from other experiments with a sufficient accuracy and 
their best values are determined from the spectroscopic 
data under question (see, e.g., [i|). Therefore, to sepa- 
rate these variables and to determine R^ one has to deal 
with at least two different transitions [l|, |^, [T2| ■ 

One of the utilized transitions is Is — 2s in hydrogen 
and deuterium, for which the distances in the interval 
between oq and 4ao are involved. As to the other transi- 
tions, a certain weighted average of data for all 2s — ns/d 
transitions is used. The value of the principal quantum 
number n for the involved excited ns/d states is sub- 
stantially larger than 2, and because of that the related 
Coulomb contribution is not larger than 1.5% and can be 
neglected as a good approximation. 

So, the Rydberg constant [1] 



R^ = 10 973 731.568 527(73) m' 



(2) 



corresponds to r ^ (1 — 4)ao. This value of i?oo can be di- 
rectly compared with another value related to a different 
distance (see below for a comparison with a value related 
to r- 10%). 

Alternatively, one can extract a value of the fine struc- 
ture constant a from the result ^ by applying the rela- 
tion 



-n-oo — 



a TTleC 
2h 



(3) 



For this purpose one should combine a value of Roo with 
h/M for a certain atom/particle and with a result or a 
chain of results, which allows one to determine M /me 
(see [l[ for detail). 

The key point in such a determination of a is a mea- 
surement of h/M. There are two high-precision inde- 
pendent experiments on that. The first measurement is 
performed on caesium atoms [l3| and the most recent 
result leads to (131 



All auxiliary data are taken from [l[. 

The other experiment has been performed on rubid- 
ium. The most recent result is [ij] 

aj^^ == 137.035 999 45(62). (5) 

The weighted average of these two results^ 

"atom = 137.035 999 59(53) (6) 

is a value of the fine structure constant related to dis- 
tances from flQ to 4ao. 



III. PRECISION MEASUREMENTS OF 

HIGHLY EXCITED STATES IN THE 
HYDROGEN ATOM (DISTANCES ~ 10^ ao) 

There is one more result on determination of the Ry- 
dberg constant, which is usually not included into the 
adjustment [l| because it is somewhat less accurate than 
the quoted above data and, what is more important, the 
result is only a preliminary one. The value 



Roo = 10 973 731.568 34(69) m- 



(7) 



is obtained [l6| from a transition between the Rydberg 
states with n ~ 30 in the hydrogen atom. Since the value 
resulted from a partial evaluation of possible systematic 
effects, to be on the safe side the uncertainty in ([7]) is 
tripled against its original value. As it was confirmed 
by the authors of the experiment [iTf , that is rather an 
overconservative estimation of the uncertainty. 

Let us now constrain an effective interaction with mass 
A, which is 

4 eV < A < 1 keV . (8) 

A contribution of an additional boson in such a case reads 



Rno — 



at r ~ ao , 
at r ^ lO^flo 



That leads to a constraint \lt 



('0.6 ± 2. 3") 



X 10 



-13 



(9) 



(10) 



That is the strongest constraint on a long-range spin- 
independent interaction, which can be derived from 
atomic physics and fundamental constants. 



ad = 137.036 0000(11). 



(4) 



^ A somewhat less accurate previous rubidium result o^j^ = 
137.035 998 83(91) from the same group [15| has not been in- 
cluded into averaging because of possible correlations. Once we 
treat that result as completely independent, the average value is 
137.035 999 39(46) (cf. lO), which does not change the 



a 



atom 

final conclusions too much. 



IV. PRECISION PHYSICS AT THE COMPTON 
WAVE LENGTH (MOMENTUM ~ niec) 

One can also compare physics at the Bohr radius and at 
the Compton wave length of an electron. A suitable tool 
for that is an examination of various results on the fine 
structure constant. We have already discussed a value re- 
lated to ao (see ([5])). Meanwhile, the most accurate value 
of the fine structure constant comes from the anomalous 
magnetic moment of an electron 



aZ}^ ^ 137.035 999 084(51) 



fll) 



obtained by combining the experimental result [19[ with 
theory ^. 

Comparing it with the value in ^ one can find a con- 
straint for 

4 keV < A < 0.5 MeV . (12) 

The results of the measurements can be presented as 



a 



a + a 
a , 



at g ^ me , 
at g ^ l/flo 



(13) 



where we refer to momentum space and assume that the 
intermediate particle is a pure vector. The constraint 
reads 111 



/= ("2.7 ± 2. g") 



a' = 2.7 ±2.9 X 10" 



(14) 



If the particle has another spin or it is a pseudovec- 
tor or a combination of vector and pseudovector (such 
as the Z boson), then the constraint is of the same or- 
der of magnitude but a factor compatible with unity can 
appear. 



V. FINAL CONSTRAINTS 

Our results are summarized in Table ID We can extend 
our analysis beyond the asymptotic limits. In particular, 
we also present in Fig. [T] res ults, which cover an interme- 
diate area of a few keV [l3 ■ 



Range 



a{r) - a{ao) 



Effect 



4 eV < A < 1 keV 



-0.6 ±2.3) X 10"" H, n = 30 



4 keV < A < 0.5 MeV (2.7 ± 2.9) x 10"" Qe - 2 

TABLE I: The constraint on the deviation of the effective 
long-range interaction a{r)/r from the Coulomb exchange 
[18| . Here, a{r) — a{oo)+a' exp{ — \r). The ao scale is related 
to hydrogen spectroscopy for transitions involving low states 
(Is, 2s). The related distance range is r = 0.5 x 10"^ m (for 



A = 4eV), 0.5 X 10" 
(for A = 0.5 MeV). 



(for A = 4 keV) and 0.4 x 10" 




A (MeV) 



lO"*^ 10"* 10"^ 0.001 0.01 0.1 



FIG. 1: Constraints on a long-term spin- independent inter- 
action from hydrogen spectroscopy and g^ — 2, including a 
constraint from hydrogen spectroscopy with low and Ryd- 
berg states (a), a comparison of low states and g^ — 2 (b). 
The lines are for \a'\ and the confidence level corresponds to 
one standard deviation. 



To extend the constraint to a region A ^ 1/ao and 
A ^ Trie, we have to modify observational equations (|9]) 
and pi|. 

Instead of ^ we find 



Roo - 

where 
-^12(2;) 



(Q!+Q'J^i2(A/(Qmc))j me 
(q+o J^3o(A/(aT?Xe))j me 



at r '-^ ao , 
at r ^ lO'^flQ 



(15) 



1 



l + x 



-2 



l + x 



2 Vl + ic 



,2 + x 
and for J^3o(A/(Q?7ie)) we have to consider a function 

2n+l / _ \ 2n-l 



(16) 



-7^30(2:) = 



2 \ _ 1 f 2 A 

2+{n+l)xJ 'r^\2+nxj 



l/(n+l)2-l/n2 



(17) 



related to n ~ 30. The transitions studied in J16l] were 
transitions between the circular states {I — n — I): 27 -> 
28 and 29 -> 30. 

The related constraint is 



(0.6 ±2.3) X 10 



-13 



-7^12 (A/ (ame)) - J^3o{X/{ame)) 
Eq. (IT^ is to be substituted for 

_ j a + a'J-g-2{X/me) , at g ~ rrie , 
1 a + a' J^i2 (A/ (ame)) , at q ^ 1/ao , 



(18) 



(19) 



where the function [2ll 



VI. SUMMARY 



J'g-2{x) = 2 -{x^ -Ax^ + 2a 



tan 



-1 /4-a;2 



^/A^^ 



2x )lnx — X 



is the function which is appUed when one calculates the 
hadronic vacuum polarization through integration over a 
dispersion variable A = ^/s. 
The related constraint is 



(2.7 ±2.9) X 10- 



Fg-2{\/me) - Ti2{X/{ame)) 



(20) 



We summarize behavior of the profile functions T in 
Fig. [5] Each of them is equal to unity for A — ;> and 
equal to zero for A — >■ cxd. The unity region is the one, 
where the Yukawa potential is indistinguishable from the 
Coulomb potential, and the effective Coulomb coupling 
constant there is a + a'. In the zero region the Yukawa 
potential shrinks to a zero-radius potential and does not 
affect any measurable values. 




A (MeV) 



FIG. 2: The profile functions, which determine coverage of the 
constraints on a long-term spin-independent interaction from 
hydrogen spectroscopy and the anomalous magnetic moment 
of an electron. The profile function J-'^o (x) is with n — 27 and 
n = 29. 

What is different in the profile functions is the charac- 
teristic value of A, which separates the unity region from 
the zero one and the shape of the transition area around 
it. 

We present in Fig. [2] four rather than three profile func- 
tions. Function J^^q is presented in two versions. That 
is due to the fact that the result Q is a preliminary one 
and the evaluation has not yet been completed. The data 
included two transitions, namely 27 — ?► 28 and 29 ^- 30 
and it is in part uncertain what their relative weights 
should be in averaging. To clarify the issue, we present 
in Fig.[5]two versions of the profile functions Tsq, namely 
for n = 27 and n = 29. One sees that the results are 
nearly identical in the logarithmic scale. Due to that, for 
the final constraint in Fig. [1] we have applied a profile 
with n = 28. 



The interaction with an intermediate particle is not 
necessarily universal. The hydrogen and deuterium data 
allow one to conclude about Uep and aed, while the 
anomalous magnetic moment deals with aee- (The un- 
certainty of the anchor data involving hydrogen and deu- 
terium is small and the accuracy of the constraint is de- 
termined in both cases above by the other pieces of data. 
That allows to reduce the ao related data to pure hydro- 
gen data.) A possible discrepancy between these two val- 
ues could be interpreted not only as related to the certain 
distance scale, but also to non-universality of one-particle 
exchange within a broader range. 

We remind that the neutrality of a hydrogen atom or a 
neutron (see, e.g., Q) is a property related to a macro- 
scopic distance and, in principle, a direct test whether 
there is an ultraweak long-range component of interac- 
tion of a hydrogen atom or a neutron at microscopic dis- 
tances is required. We will address this problem else- 
where. 

For the b constraint in Fig. [1] we have utilized two 
values of the fine structure constant, namely those from 
(7e — 2 and from atomic physics. In principle, one can use 
other less accurate values (see Fig. [3] and [l| for a review 
of data and references) to compare with ag-2 ■ That pro- 
duces a constraint similar to the b line in Fig. [TJ It is 
weaker for A '3> ciq^ , but has the same profile behavior, 
while for A < Oq the limitation increases slower than 
the b line. However, as we see in Fig. [U in this area the 
a constraint is already applicable and such extension of 
the b constraints to the lower mass range cannot improve 
the whole constraint. 

Another issue is model dependence. The substitute ([T]) 
is to constrain a situation with a single light intermedi- 
ate particle. If there are a few of them and their cou- 
pling constants are comparable, the situation becomes 
more complicated. Still, if the related coupling constants 
are of the same sign, the results of the paper are still 
applicable. However, if the Yukawa terms involve inter- 
actions with different radii, while the coupling constants 
have compatible absolute values and opposite signs, the 
constraint derived is completely out of validity. As an 
example, one can consider a long-range correction of the 
form 



^Air 



3-Air 



where the masses are of the same order (e.g., A2 = 2Ai). 
Indeed, in this case the sensitivity of atomic energy levels 
to a' is substantially reduced in most of the area, where 
the constraint in Fig. [T]is effective. 

Concluding we emphasize that atomic physics allows 
one to access distances much larger than the Bohr ra- 
dius and our limitation on a' is valid for distances up 
to 0.02 /im. That can be compared with the limitation 
on the fifth force from Casimir-effect experiments |24 |. 
which also reached the range below 1 /zm. Comparing 



R^(NPL-f 



r_(lo)(NIST-89) 



r' „(lo)(NIM-95) 



R^(NIM-95) 



1— •— I R^(NIST-97) 

1 R^(NMI-97) 

j„(lo)(KR/VN-98) 

— • 1 R^(LNE-01) 

m h/m (Cs) 



■ h/m (Rb) 
• a^(UWash-87/Th06) 
• a_^(HarvU-06/Th06) 



FIG. 3: Determination of the fine structure constant a in the 
CODATA-2006 adjustment [3|, where the references can be 
found. The vertical band stands for the adjusted CODATA- 
2006 value. The figure is reproduced from [23] . 



it, we find that our constraints starting from a distance 
below 30 nm, where the Casimir-effect constraints stop, 
are weaker by a few orders of magnitude. However, they 
are extended to shorter distances, not accessible for the 
Casimir-force related experiments. 

What is more important, our constraint is complemen- 
tary to the Casimir-effect limitation. The latter deals 
with the neutral bulk matter only and cannot constrain 
a massive ultraweak photon 7*. In contrast to that, our 
constraint in Table |T] and in Fig. [T] is covering such a 
case. Our results for the coupling constant a' related to 
7* are consistent with zero and the limitation at atomic 
distances is in the range of a few parts in 10^'^ to a few 
parts in 10^^ depending on details as discussed above. 

The keV mass range has been partly explored by means 
of astrophysics and cosmology. Such constraints [a, |23 
mostly deal with real particles and their propagation 
through the matter. Details, such as the lifetime, various 
couplings and collision rates, are involved and they may 
be related to results in terms of a' and A only through 
certain mo del- dependent relationships. Our constraints, 
free of such model-dependent suggestions, are comple- 
mentary to those. 

Another option, also not covered by Casimir effect is 
due to a spin-dependent long-range interaction to be con- 
sidered elsewhere. 



the strength, we note that in the case of Casimir ef- 
fect the conventional parameterization for the Yukawa 
potential is Gmira2ae~^/^ /r^ where rrii is a mass of 
a bulk body and G is Newtonian gravitation constant. 
There is no direct mo del- independent correspondence 
between mass-related and charge-related parameteriza- 
tion. The most important part of the conversion fac- 



tor is Gr. 



6 X 10 (since in bulk matter there is 



roughly 0.5-1 charge, e.g., baryon or lepton charge, per 
nucleon-the rough factor of 0.5 comes from the fact that 
we have an electron per from 2 to 2.5 nucleons). Using 
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